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Abstract 



We write the probabilities for neutrino oscillations in uniform-density matter 
exactly in terms of convention-independent vacuum neutrino oscillation param- 
eters and the matter density. This extends earlier results formulating neutrino 
oscillations in terms of matter-, phase-, and trace-invariant quantities. 
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1 Introduction 

The advent of oscillation experiments in which neutrinos traverse significant distances 
in terrestrial matter before observation make the relationship between the vacuum 
oscillation parameters and the oscillation observables in matter of considerable inter- 
est. Some of us recently pointed out PP0 the significance here of matter-invariants: 
neutrino oscillation parameters which are invariant under the influence of matter may 
be used to simplify the relationship between neutrino oscillation observables and the 
vacuum neutrino oscillation parameters. 

The effective neutrino oscillation Hamiltonian in the flavour basis (the weak basis 
in which the charged lepton mass matrix is diagonal) is given by H = MM' /2E, 
where M is the neutrino mass matrix and E is the neutrino energy. In vacuum, this is 
diagonalised by the conventional MNS mixing matrix [3], U: U^HU = diag(Ai, A 2 , A 3 ), 
where the Aj = mf/2E are the vacuum eigenvalues. The effects of matter on the 
propagation of neutrinos is described by the addition of the Wolfenstein term, a = 
\f2GpN e (N e is the number density of electrons in the matter), to the H ee element: 

H = H + diag(a, 0, 0) (1) 

which modifies the eigenvalues and the elements of the MNS matrix in a non-trivial 
way. We denote matter-modified parameters by quantities with a ~. 

The electron density in matter is not C-, CP- or CPT-symmetric: for antineutrino 
propagation, a in Eq. keeps its magnitude but changes sign. Thus, the matter- 
modified eigenvalues and mixing matrix for antineutrinos are different from those for 
neutrinos. We treat the general case, leaving the matter density, a, a free parameter, 
and comment further on the relationship between the neutrino and antineutrino cases 
in Appendix A. 

The formula for the appearance and survival probabilities as a function of propaga- 
tion distance, L, when neutrinos pass through uniform density matter may be written 
in its usual form, but in terms of the matter-modified parameters as follows: 

P{v a ^vp) = 5 a/3 -4^^sin 2 (A^L/2) 

+ 8 J aP sin (A 13 L/2) sin (A 23 L/2) sin (A 31 L/2) (2) 

where the 

K% = ReiU^U:^), (3) 
parameterise the magnitudes of the T-even oscillations and 

J a p = ImiUaiU^Upj) = ±J (for a ^ (3) or =0 (for a = (3) (4) 

parameterises the magnitude of the T-odd oscillations 1 . The eigenvalue differences in 
matter, Ajj = (rhf — m^)/2E may be calculated in terms of the vacuum parameters, 

1 We prefer the "T-even" and "T-odd" labels to the "CP-even" and "CP-odd" ones, since matter 
introduces an extrinsic CP-odd contribution into the intrinsically CP-even terms, while spherically- 
symmetric matter profiles respect T-invariance for neutrino propagation between points at equal radii 
(eg. on the surface of the Earth) . 
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Ajj and U a i, and the Wolfenstein term, a, using the solutions of the cubic characteristic 
equation of the matter-modified Hamiltonian jl] [3] . The matter-modified MNS matrix 
elements, U a i, may be similarly calculated, but are rather complicated functions [3] of 
the vacuum parameters and the Wolfenstein term. It is the aim of this paper to simplify 
as much as possible the relationship between the observable oscillation amplitudes in 
matter, ^K % h and 8J, and the vacuum parameters. 



2 Matter Invariance 

The idea of matter-invariance is based on the observation that all quantities H a p 
in Eq. (0) other than H ee are matter-invariant, and appropriately combined, can be 
related to observable parameters (H ee , because of its trivial transformation in matter, 
may be said to be "matter-covariant"). The first application of these ideas showed that 
Jarlskog's determinant jH] was matter-invariant and led to the so-called NHS relation 

m 

A12A23A31J = Im{H^H^ T H re ) = lm{H etl H^ T H re ) = A 12 A 23 A 31 J. (5) 

This relation was used to write exactly the T-violating part of the matter-modified 
oscillation probability, Eq. (j2J), in a very simple and compact form: 

P 7 iy a - v p ) = 8 t 12 t 23 t 31 J sin ( A i2 L /2) sin (A 23 L/2) sin (A 31 L/2) (a^f3). (6) 
Ai2A 2 3A 31 

The advantage of this formulation is that it does not require the matter-modified 
MNS matrix elements, whose expressions in terms of vacuum parameters are quite 
complicated |5J- The matter-dependence is confined to the eigenvalue differences, Ay, 
for which it is somewhat more straightforward 

The matter-invariant approach was subsequently extended [2] to the T-even part of 
the oscillation probability, for which analogous (but less simple) expressions to Eq. (jSj) 
were found. Analogues of Eq. (jUJ), valid in approximation were obtained, but exact 
formulations were not. 

The matter-invariant approach was next applied ingeniously [7j to provide, after a 
lengthy derivation, exact expressions for the T-conserving coefficients [S] in terms of 
the effective Hamiltonian elements and the matter-dependent eigenvalues: 

T^j \H a p\ 2 \i\j + \Q afS \ 2 + Re(H a f 3 Q* a/3 ){\i + \j) 

K «P= A A A A {a ^ P) - ( ) 

Here Q* af3 is the cofactor of H a p and is matter-invariant for a or f3 = e, while H a p is 
matter- invariant for a or f3 7^ e. For the two off-diagonal cases with {a, [3} = {e, /i} or 
{e, r} (the K^n are symmetric under the interchange of a and ft), this approach isolates 
the matter-dependence to the A, alone, in the same spirit as, eg. Eq. (JHJ isolates it to 
the A^. However, for the case {a, (3} = {fi, r}, Q a @ is not matter-invariant, which 
complicates the formulation, and spoils the symmetry of the formulation between the 
different components. 
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Neutrino oscillation observables can depend only on the differences between mass- 
squared eigenvalues and so must be "trace- invariant," ie. invariant under transforma- 
tions in which any multiple of the identity is added to the effective neutrino Hamilto- 
nian; a change in the trace is equivalent to a change in the overall phase of the neutrino 
propagation amplitude. Observables are, in addition, "phase-invariant", ie. invariant 
under phase transformations of the neutrino mass eigenstates. Hence, it must be pos- 
sible to write the relationship between the observables and the vacuum parameters 
entirely in terms of trace- and phase-invariant quantities. Eqs. (JHJ) and (JHJ) are ex- 
amples of this. While the particular combination given in Eq. ((7J) is of course both 
trace- and phase- invariant, this formulation suffers the difficulty that neither the Q a p 
nor the \ are trace- invariant. These individual quantities cannot therefore be related 
to observables of neutrino oscillations. Before the values of the Q a p or the \ can be 
specified, an artificial offset of the neutrino masses must be chosen. In the applications 
cited in [7] the offset is arbitrarily set so that m\ = in vacuum. 

In the remainder of this paper, we provide a unified formulation, using matter- 
invariants which are trace- and phase-invariant. The matter-dependence is isolated 
in factors which depend only on the eigenvalue differences, A^-, and the matter den- 
sity itself. We find the exact T-even analogues of Eqs. (J3J) and ©, and hence exact 
convention- independent, matter-covariant expressions for the observable neutrino os- 
cillation probabilities in terms of vacuum parameters. 



3 Matter-Covariant Derivation of Oscillation Probabilities in Uniform Den- 
sity Matter 

We provide a matter-covariant derivation of the neutrino oscillation probabilities given 
in Eq. (J2J) . The amplitude A a/ 3 for a neutrino of flavour a to be detected as a neutrino 
of flavour (3 in matter of uniform density is given as a function of propagation distance 
L by the (matrix) equation: 

A = exp(-iHL) (8) 

where H is the effective neutrino oscillation Hamiltonian of Eq. (JTJ. The general 
theory for a function of an operator jHj enables the exponentiation to be performed 
directly in the flavour basis: 

A = J2X i exp(-t~X t L) (9) 

i 

where the Aj are the eigenvalues of H and the Hermitian projection operators X 1 are 
given by: 

X 1 = FWff-A,-) (1Q) 



(j^k^i). (11) 



The second equality is specific to three families of neutrinos. Comparison with the more 
conventional approach in which the Hamiltonian is diagonalised before exponentiation 
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allows us to identify the elements of the X % with the familiar combinations of the 
lepton mixing matrix elements |lUj : 

Xap = UaiUpi (12) 

(no summation over i is implied). The X^o for a ^ (3 are not phase- invariant, and 
are therefore not observables (although the diagonal elements X l aa = \U a i\ 2 are). All 
components of the X % are however trace-invariant, this property being manifest in the 
combinations {H — A*) which appear in Eqs. (fTUj) and (fTTj) (or in the fact that diago- 
nalization of H, yielding the U's appearing in Eq. ([12)1. is a trace- invariant process). 

The neutrino oscillation probabilities of Eq. (J2J) are given by the squared amplitude 
P{y a — > vp) = \A a /3\ 2 , which contains real and imaginary projections, K^g and J 
respectively, of the phase-, trace-, and convention-independent plaquettes fTTj : 

u a v;Pifi Pj = K P x j : P (i^j). (is) 

Comparing Eqs. (fT2"j) and (JT3j) . one sees that the X l aj3 represent an intermediate calcu- 
lational step between the U a i and the observable plaquettes. Unlike the U a i and the 
plaquettes however, the X^ have a simple unitarity relation: J2i^a/3 = <W- The X's 
close the relation among the fundamental parameters H and A: the relations 

\ = Tr(X l H) , (14) 

and 

(H-~X k ) = J2^ikX i for given k, (15) 

i 

together with Eq. (fTTj) demonstrate the equal status of the X, H, and A. Any pair of 
these sets of quantities encapsulates equivalent information. 

We now proceed to develop Eq. (fTTj) by explicitly calculating the X l af3 in terms of 
the H and A elements. Following this, we relate the matter and vacuum values of the 
Xs. From Eq. (fTTj) . we have for the case a (3, 

AjjAjfcX^ = [H 2 — (Xj + \k)H] a /3 

= Ha-yH^ - (# 77 - \i)H a/3 (a ^ (3 ^ 7, i ^ j ^ k) (16) 

where we have used the fact that Tr(H) (= T) = J2 a H aa = Eq- (fTUj) was 

obtained in [7j, in a less straightforward manner. We note that on the right-hand side, 
all factors are matter-independent except for (if 77 — Aj). This factor is also problematic, 
as it contains the Aj, which are not directly observable in neutrino oscillations. We 
deal with the both problems by the substitution in terms of observable quantities: 

(# 77 - A;) = (# 77 - ~T) - Aj + cw 7 (17) 

where the 

A i = -(A i ,- + A ii; ) (j^k^i), (18) 
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are the eigenvalues of the reduced matter-dependent Hamiltonian (H — |T), and t> 7 
is the "7" element of the vector: 

*=<!-4-5> (l9) 

which breaks the symmetry between the three lepton flavours for a 7^ 0. The matter- 
dependence in Eq. (IT7jl has been isolated to its last two terms. Substitution from 
Eq. (fTTjl in Eq. (fTHjl shows that the quantities 

R a p = AijA ik X % a(3 + (av 1 - Ai)H a p, {a ^ (3 ^ 7, i ^ j ^ k) 
= — (H 77 — -T)H a p 

= H a7 H 7 p — (iT 77 — -T)H a p (20) 

are matter-invariant, ie. R a p = R a p for all a 7^ (3. 

We can now relate the matter values X 1 and the vacuum values X 1 : 

-^ a/3 = A A ( a f ^ f Tl J / M ! )- ( 21 ) 

The matter-invariants appearing in Eqs. (j20j) and (J2TJ) may themselves be expanded 
in terms of vacuum values of As and Xs: 

-. cyclic k fixed 

H a p = H a p = - £ (Ajj + A ifc )X^ = £ A lfc A^ (a ^ 0) (22) 

j cyclic -j A; fixed 

Rap = Rap = ~ E AyA^ = 3 E ^ fc (Ay + Ay)^ (« ^) (23) 

where the former is obtained by summing Eq. (|15|) over (with a = 0) and using 
unitarity, and the latter by summing the first line of Eq. (J2*U|) over % (again with a = 0) 
and using the fact that T,i^i = 0. Using Eqs. flU or (j2U) with Eqs. (221 and (J2BD 
allows the matter- modified values X l af5 = 11^11% to be calculated without solving for 
the U ai themselves and without any dependence on the offset of the eigenvalues. 

We remark that the X^p, X^p, H a p and R a p (a ^ (3) are not observable, all having 
a common phase-convention dependence. Rather, the X l a g are to be considered as 
building blocks of the observables, K l h. The diagonal components, X' L aa = \U Q i\ 2 , are 
however observable; their calculation by a similar method is discussed in Appendix B. 

It is now easy to calculate the T-conserving and T-violating oscillation coefficients 
K^g = Re(X l a/3 Xap) and J = Im^X^X^) in similar terms and exhibit their matter- 
dependences. From Eq. ()20)1 we find: 



f>ij _ Ay \H a p\ + ff 7 Re(H a pR* af3 ) + \R a p\ 

^ap~ A A A A [a * Ph 1 ' 
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where all terms and factors are independently observable and the matter-dependence 
is confined to the coefficients 



_, i_„ i 2 

A* = A t Aj + a Vj A k + a> 7 ) 2 = g(-2Aj + A ki A kj ) + -a Vl {A kt + A kj ) + -a> 7 + -), 

(25) 

= -(A k + 2a Vl ) = --(A w + A kj ) - 2av 7 (i ^ j ^ k) (26) 

and the denominator. There is thus no dependence on the matter-modified mixing 
matrix elements, the matter-dependence entering only via the explicit a-dependent 
terms and the Ay, which are given by standard expressions in terms of vacuum pa- 
rameters and the matter density in Eq. (|24j) is similar to the exact formula for 
the T-even oscillations, Eq. (J7J), except that here, Eq. p^ty * s composed of explicitly 
observable quantities, and the matter-dependence has been isolated for all a ^ (3. 
Similarly 

- (Aj - A 3 )Im(H a pR* a p) 

J = — = — = — = — (a f: B) 

A 12 A 23 A 3 iA ij 

Im ( H efJL H^ T H Te ) 
A12A23A31 

where we have used the relations 



(27) 



H a f3R* a/3 = (H e ^H^ T H Te )^ — (if 77 — -T)\H a/3 \ 2 (28) 

and 



(A t -A 1 )=A tl . (29) 

Eq. (j2Zj) is simply the well-known result which leads to the NHS relation [1J |7| 
of Eq. © above. 

4 Exact Oscillation Probabilities in Terms of Vacuum Parameters 



Using Eq. (|24jh we can now solve for the matter-dependence of the K^ g in terms of 
A 12 A 23 A 31 A^ = A 12 A 23 A 31 A ij K% + ^ p (a^(3), (30) 



the vacuum K%g 



where 

k% = (i* - A*)\H a p\ 2 + (J3* - B>;)Re(H aP R* ap ) {a^^^i^j^k) (31) 

with A* and B* given in Eqs. (J2SJ) and (JHJ) and A k = |(-2A?. + A ik A jk ) and 
B* = — |(A fc j + Afej) (which do not depend on 7). Equation (|3Uj) is the exact T- 
even analogue of the T-odd invariance, Eq. (jSJ). It is slightly more complicated in the 
sense that the matter-modified values K% differ from the by an inhomogeneous 
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term K^/Ai 2 A23A 3 iAj ;? -, as well as by a scale factor. The K^p clearly vanish in the 
limit a — > 0, as they should. 

In order to complete our specification of the K^* in terms of the vacuum K^s, it 
is necessary to give \H a/3 \ 2 and Re(if Q/ gi?* /3 ) in Eq. (J3"T]) in these terms. The former is 
given by 0: 

cyclic cyclic 

- £ AJ>% H#*l a = M = - £ A "A (32) 

(v) (*i) 
while the latter may be derived from Eqs. (J22j) and as discussed in Appendix C: 

cyclic 

MHaflKp) = -3 E + A ^)^S (« ^ /*)■ ( 33 ) 

Eqs. f)30 |l -([3^ )l are the main result of this paper, which, along with Eqs. (|23j) and (J2HJ), 
allow the if^g to be calculated in terms of the vacuum parameters Ay and if^, the Ay 
and the matter parameter, a, thereby avoiding the need to find the matter-modified 
MNS mixing matrix. These formulae are, furthermore, all convention- independent. 

In Fig. (JTJ) we plot the nine K^n, a 7^ /3, i 7^ j, for neutrinos and antineutrinos 
traversing the Earth's mantle, as functions of the neutrino energy, calculated using 
Eqs. flU-flSBI). We take Am 2 12 = 5.0xlO" 5 eV 2 , Am 2 3 = 2.5xKT 3 eV 2 , sin#i 2 = 0.58, 
sin 823 = 0.71, sin ^13 = 0.05 and 5 = tt/4. The similarity of the first and second rows 
reflects the approximate v v -v T symmetry [13] of the vacuum MNS matrix. Except 
near a matter-resonance, the smallness of the values in the upper-right 2x2 sub-block 
reflects the smallness of |£/ e 3|. 

From Eq. (J2J), we can write the exact expression for the appearance probabilities 
in uniform density matter: 

P{v a — ► vp) = \A a p\ 2 = 

+ 

(a ^ where Ky is given by Eq. (|3T|) in terms of vacuum quantities and the and 
B*, which in turn depend only on a and the Ay. Although we have calculated only 
appearance probabilities in Eq. (|34p. survival probabilities are calculable in similar 
terms directly from them using unitarity. This completes our derivation of the exact 
matter-covariant formulation of neutrino oscillation probabilities in uniform density 
matter, in terms of vacuum oscillation parameters and the matter density. 

Our formulae also hold for antineutrino oscillations. For antineutrinos, the signs of 
J and of a are opposite to those for neutrinos. These sign changes alter the effective 
Hamiltonian, and the values of the eigenvalue-differences, Ay, are changed in our 
formulae. In Fig. ^ the consequent differences between the antineutrino and neutrino 
oscillation coefficients are clearly seen. We discuss this more in Appendix A. 



-4E 



a/3 



A12A23 A31 Ay K%g + K 3 

Ai2A 23 A 31 



A 12 A 2 ,A,iA, : 



sin 2 (AijL/2) 



A12A23A31 



12^*23 ^31 ^-Hj 

J sin (AiaL/2) sin (A 23 L/2) sin (A 31 L/2) 



(34) 
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Figure 1: The nine coefficients K^g, a ^ (3, i ^ j, for neutrinos (dashed lines) and 
antineutrinos (dotted lines) traversing the Earth's mantle, as functions of the neutrino 
energy. We take Amj 2 = 5.0 x KT 5 eV 2 , Am 2 3 = 2.5 x 1(T 3 eV 2 , sin6»i 2 = 0.58, 
sin 9 2 3 = 0.71, sin# 13 = 0.05 and 5 = 7r/4. The corresponding vacuum values are 
indicated by an arrow. 



Appendix A: Matter-Effects for Antineutrino vs. Neutrino 

Since matter is inherently CP-violating, it affects antineutrinos differently than neu- 
trinos. This can be summarised by noting that the sign of a is negative for antineutri- 
nos, positive for neutrinos. Thus, any matter invariant is also a neutrino-antineutrino 
invariant, while any matter dependence breaks the neutrino-antineutrino symmetry. 
The explicit breaking of this symmetry is readily obtained from some of our derived 
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formulae. From Eq. (J2Tj) we get 



([Ai] v - av 1 - [Ajjp - av^)H a/3 (a ^ /3, j ^ k ^ i) 



-([A^ + A^]^ - [Ajj + A ifc ] p ) - 2<w 7 [> if Q/3 . 



From Eq. (J27j) we get 



A 12 A 23 A 31 J 



A 12 A 23 A 31 J 



Similar but longer expressions may be written down for 



Ai 2 A 23 A 3 iA ij ii'^ g 



A 12 A 23 A 31 A ij K^ p 



K 



a/3 



K 



a/3 



using Eqs. (f3T)J) and (f3~TJ) (the terms quadratic in a cancel) and for 



Aij Ajfc | L^jj | 



AijAik\U a 



using Eq. flUJ. 



(35) 



(36) 



(37) 



(38) 



Appendix B: Matter- Vacuum Relation for |E/o.i| 2 

The starting point is the expansion of the diagonal element of Eq. (fTTj) : 

~A t3 A lk T aa = ]T (H — Xj) a(T (H - Aj;) CTQ (j^k^i) (39) 

= (^-A.O^-A^ + I^^ + I^I 2 (40) 
Inputting Eq. (fTTj) for (H aa — Aj) then isolates the matter dependence. The result is 



Ajj Ajfc | U a i | 



A^A^U^ 2 + (il - A l a ) - (B« - B l a ){H aa - -T) 



(41) 



where the coefficients are the same as those in K l £p, Eq. (jHIJ. The factor {H aa - §T) 
can also be put in terms of our vacuum observables using the appropriate diagonal 
component of the matrix equation: 



1 



cyclic 



(H--T) = -Y: (Ay + A lk )X* 
which is obtained by summing Eq. (fTo'j) over k (taking the vacuum case). 



(42) 



Appendix C: Matter Invariants in Terms of K^p and |t/ a i| 2 

The derivation of Eq. ()33)1 follows straightforwardly from Eqs. (j22J) arid (|2*3*j) . utilising 
the useful relations between the and the |X^| 2 : 



I V* 1 2 



(43) 



10 



and 

K% = \{\X k a ^-\X^-\X^) Va^, (i^j^k) (44) 

which are themselves easily derived from the unitary condition Y^iX^* = 0. 
These can also be used to find the matter-invariant 

-i cyclic 

l^| 2 = -o E A 2 .(A ifc + A ife ) 2 i^ (45) 

which, in addition to Eqs. (J52j) and (jHHj) . completes the set of matter- invariants used 
in Eq. (J23j) ( we did not need to specify these in the main text, because we used the 
substitution of K l K instead to find Eqs. (J3UJ) and (pTTjO . 

We can also use Eqs. (J43j) and (J44j) to find the set of three matter invariants used 
in Eq. instead in terms of \X l a/3 \ 2 = \U ai \ 2 \Upi\ 2 : 

cyclic 

\H a p\ 2 = ]T A^A lk \X l J 2 (46) 

-i cyclic 

Re(H a pR* ap ) = - £ AyAttCAtt + A, fc )|X^| 2 (47) 

2 cyclic 

\R a p\ 2 = - £ A ii A ife (A ii A ift - 2A 2 fc )|X^| 2 . (48) 

Substituting these into Eq. (J23j) yields a formula for the which depends only on 
the moduli-squared, |£7 a j| 2 , of elements of the vacuum MNS matrix (in addition to the 
Ay and the matter-density, a). This alternative formulation may be considered more 
convenient to use than the one using the vacuum K 1 ^. 
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